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(57) Abstract: A method of obtaining data for use by a receiver of a satellite positioning system or a GNSS comprises deriving the 
data remotely &om the receiver by a server (200), using meteorological information and a regional or global three dimensional map 
of grid points from which it computes tropospherical delays by ray tracing through the rcfractivity field derived from atmospheric 
measurements of pressure, temperature and water data content, such measurements being available from meteorological bodies. 
When used to enhance position determined by a user receiver that includes a non-meteorological, climate based model (130) giving 
zenith delays and means (130*) to map them to particular inclinations, the server also includes a copy of such non-meteorological 
model (230) and provides its ray traced delay values as zenith delays. The sets of zenith delay values for corresponding grid points 
are compared in the server (260) and modifications developed (preferably in fractional form) by which the non-meteorological delay 
values require correcting to be accurate. The correction sets are reduced by image compression techniques (270) and transtnitted via 
the satellites (1 lOi etc) of the GNSS at low data rate to the user receiver, which receiver simply applies the corrections to the Zenith 
delays derived by its own model- If a user position is known, the server may derive accurate tropospheric delay values dircctiy for 
the receiver position direcdy for transmission. 
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CORR£CTION OF TROPOSPHERE INDUCED ERRORS IN GLOBAL POSITIONING SYSTEMS 

The present inv^tion relates to developments intended to increase the accuracy 
obtainable from global navigation satellite systems (GNSS). 

5 At present there are two publicly available GPS systems, known as NAVSTAR, 
owned by the USA, and GLONASS owned by the Russian Fed^ation. These have 
been in existence for arouid two decades, but in the near future it is hoped that the 
E\iropean regional augmentation of GPS will start to provide its services, followed 
within a few years by a European system under the name of GAULEO. 

1 0 The existing systems have been progressively refined so that using a differential phase 
implementation a locational accuracy of less than 2 cm can potentially be achieved 
over a baseline of 1000 km, but with a cost in computation and in tiie time taken to 
determine the location. Real time or near real time measurements have a 
correspondingly lower resolution, and at present the requirements for high precision 

15 mean that additional augmentations are necessarily employed to supplement the 
GNSS information. Furthermore, these could include a receiver taking measurements 
from many satellites, up to all those visible to it whereby to calculate an over- 
determined position solution and rejecting inconsistent data to improve the accuracy 
of the position solution. Such a system may use data from more than one constellation 

20 of GNSS satellites, GPS and GLONASS. 

Although GNSS is used mainly for establishing the location of a user having a 
suitable, usually mobile, receiver, it is also used in respect of providing accurate time 
signals to users whose locations are already known or do not need to know. Singile 
user position determining sets have simple receivers of satellite transmissions and 
25 circuitry that effects modelling of at least some atmospheric effects that influence 
signal reception so as to go some way towards eliminating errors in calculated 
position. 
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However, whether fhe user is interested in obtaining a position or a time 
measurement, a significant error arises from the inability to model accurately the 
delay to the GNSS signals caused by the atmosphere, namely the ionosphere and 
troposphere. 

5 Satellite navigation users generate their three-dimensional position and time solution 
by processing four 4 or more pseudorange measurement to four or more satellites, A 
pseudorange measurement is the difference between the satellite clock time at signal 
broadcast and user receiver clock time at reception. The pseudorange observation is 
therefore related to the radio propagation time and therefore range between satellite 
10 and user. As estimates of the satellite position are known (they are broadcast by the 
satellite) a user can solve for the four unknowns (three-dimensional position and 
time) using four or more pseudorange observations. As part of the user's 
navigation/time solution filter pseudorange observations are corrected for variations 
in radio propagation time from that of fiee-space propagation. 

15 

In the user's navigation/time solution filter^ a number of corrections are applied to the 
raw pseudorange measurement including tropospheric, ionospheric and relativistic 
corrections. 

20 It has been suggested in WO-Al-03/069366 how to accommodate ionospheric delays 
and by use of a so-called server site which receives GNSS satellite signals, derives 
correction factors applicable to GPS receivers in tiie vicinity before broadcasting them 
locally so as to be received by such a GPS receiver and used to modify the on-board 
model used to correct such delays. For ionospheric delays, which comprise a small 

25 degree of signal path refraction and a more significant change in signal velocity, the 
dela>^ and corrections th^efor are substantially constant over a period of time that 
requires updating of correction data at most a few times per day. 

Tropospheric ejSects on the other hand are relatively fast changing (or short-lived) and 
geogr£q)hically localised, resulting primarily &om weather or meteorological 
30 phenomena rattier than climatic phenomena. However, the troposphere constitutes 



wo 2004/095055 PCT/GB2004/001676 

3 

one of the largest identified sources of error in ttie effect that it has on signals 
propagating therethrough. The troposphere introduces ray bending and therefore an 
increase in signal path that constitutes a signal delay which is influenced by a number 
of meteorological factors, but particularly water content. Tropospheric delays are 
5 dif&cult to model shnply. 

Traditionally, tropospheric delay has been handled by the use of global tropospheric 
delay models that work firom so-called climate parameters that relatively invariant and 
can be stored in the user receiver, but these parameters at best constitute an average or 
seasonal expectation, but not one that is meteorologically based, that is, based upon 
1 0 current, recent or predicted weather conditions. 

One such model that is used and may be built into a portable GPS receiver is the 
RTCA tropospheric zenith delay model for WAAS users described in "Minimum 
Operational Performance Standards for Global Positioning Systems / Wide Area 
Augmentation System Airborne Equipment" RTCA D0229C, November 200 L 

Such model is useful insofar as it simplifies tropospheric delays to zenith values 
(identified herein as or ZTD) but there is still the need to map these for elevation 
effects caused by low satellite inclinations to the user. One such mapping model is 
described by Niell in "Global mapping functions for the atmosphere delay at radio 
wavelengths" Journal of Geophysical Research Vol 101, No B2, Pages 3227-3246, 
February 1996. 

However, although these models permit incorporation into a user receiver they are 
inherently limited in ability to accommodate changes in tropospheric conditions that 
affect signal delays caused by the constantly changing, and localised weaflier. 

Although models exist for deriving accurate tropospheric data by taking into account 
2 5 the meteorological conditions in one or more regions, such as by numerical weather 
prediction (NWP), the localised nature and thus large amount of data generated has 
been perceived as confirming that presently they cannot be used to sensibly improve 
upon practicable devices; that is, due both to this data being too large to be sent over 



15 



20 



wo 2004/095055 PCT/GB2004/001676 

4 

communication systems that are available to mobile users and the limited capacity for 
processing within a reasonable amount of time. 

The present invotition provides a method of obtaining tropospheric delay data for use 
in a satellite positioning system or GNSS comprising the steps of generating for a user 
5 location, at a location remote jfrom the user location and from meteorological 
information, at least one accurate tropospheric delay value, applicable to the user 
location for comniunication as a troposph^c delay correction to a said user. 

Preferably said accurate tropospheric delay values are derived by a ray tracing 
technique. The accurate tropospheric delay values may be derived by three- 
10 dimensional refractive index field generation. Furthermore, it is also preferred that 
said meteorological information is based on numerical weather prediction (NWP) 
data. The meteorological model or each said tropospheric delay value correction 
derived therefor may be augmented by directly observed meteorological data. 

In an embodiment of the invention, appUcable to a user whose position is not 
15 accurately known, the method may comprise generating, from a first model which is 
known per se, a first set of approximate tropospheric delay values applicable to 
various user geographical locations, generating from a meteorological model 
employing such meteorological information a second set of tropospheric delay values 
that are accurate and applicable to said various user geographical locations, 
20 developing a set of delay value modifications for use with said first model so that it 
can provide a set of tropospheric delay values substantially in agreement with &e 
second set, and expressing the set of modifications as a set of tropospheric delay 
corrections for communication to a said user. 

The first model is based on non-meteorological parameters, which parameters 
25 comprise at least one of time of year, latitude and altitude. The non-meteorological 
parameters may fiuther comprise at least one of longitude and time of day. 

In the method of this embodiment the first and meteorological models develop sets of 
tropospheric delay values comprising zenith tropospheric delays. The first model may 
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contain a mapping function relating tropospheric delay at a given elevation angle to 
the zenith tropospheric delay, and said set of delay value modifications may comprise 
a set of modifications for use with the mapping function of the first model. 

Preferably, the modifications to the delay values are the differences between 
5 corresponding values of the sets attributable to the first and meteorological models. 
The corrections to be communicated may be the modifications per se or, preferably, 
the modifications expressed as a fractional change from the values of the first set, for 
example as a percentage. 

Thus a correction may be effected as an addition or multiplier to any value generated 
10 by use of the first model. 

The accurate tropospheric delay values are derived by a ray tracing technique, to 
determine the path of the satellite signals through the troposphere to tiie user and 
hence estimate the delay from a direct path, and possibly employing three- 
dimensional refractive index field generation. 

15 It is possible to effect comparable corrections to mapping functions in such a first 
model that also rely upon paths affected by tropospheric delay. 

The meteorological model may be based on numerical weather prediction (NWP) data 
for a region of the earth or real time meteorological data or both. In particular, the 
meteorological model or each said tropospheric delay value correction may be 
20 augmented by directly observed meteorological data, such as available in two- 
dimensional form from some imagiug satellites. 

The region for which data values are obtained may be substantially global or may be a 
smaller region as defined in NWP schemes as mesoscale maps. 

In both cases it is possible to generate a set of (zenith) tropospheric delay values for 
2 5 each of a grid of locations over a region, as a two-dimensional array defining the 
geographical points at each of which is a delay value. Thus it is possible to create 
from the meteorological model and non-meteorological model zenith tropospheric 
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delay value modifications and a set of corrections as a data array having values 
determined for individual grid points on the earth's surface, and the set of values 
comprises a distribution of said modifications over at least part of the earth's surface. 

Having regard to the nature of the delay values and corrections, data reduction may be 
5 applied to the correction set, deriving a reduced data set for cocomunication to a user. 

Preferably, insofar as the nature of the correction set produces a data file analogous to 
that of a greyscale image reduction of the data size of the correction set is 
accompUshed by an image compression process, conveniently, but not essentially, by 
lossy data reduction such as according to a wavelet-based JPEG 2000 or cosine JPEG 
10 standard. 

By effectmg a suitable level of data reduction, it becomes possible to communicate 
the correction set data to a user over a communication or data channel of limited 
bandwidth. 

Thus haviag created a correction set suitable for use by a remote user having a 
15 satellite signal receiver, in accordance with this embodiment of the invention at least 
part of the correction set may be communicated via at least one orbiting satellite by 
transmitting the correction set as a reduced data image file to a said satellite and re- 
transmitting at least part of the set to a user firom a said orbiting satellite. 

Preferably this is achieved by communicating the corrections to at least one orbiting 
2 0 GNSS sateUite from which user receives signals to establish at least one of position 
and time. In order to reduce even fiirther the amount of data to be transmitted or re- 
transmitted, the method envisages communicating the image data to a said satellite for 
re-transmission of only that part of the correction data that can be of use to a user in a 
region within range of said satellite. This may be achieved by transmitting only said 
2 5 part of the correction data to the satellite or transmitting all of the data but causing the 
satellite to re-transmit only said part. 

hisofar as such a satellite has limited capacity to transmit signals additional to those 
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akeady transmitted and is in general only able to transmit any data, including 
correction data, periodically, it must be borne in mind that the meteorological 
environment is changing continuously as weather features vary their position in 
relation to the mapped region. Thus in addition to deriving tropospheric delay values 
5 associated with grid points of the mapped region it is necessary to apply data 
reduction sufficient to permit transmission of all or part of a said corrections useable 
by a user within a time, dictated by transmission availability and transmission rate of 
the satellite, substantially lower than the validity time of the meteorological 
information used by the meteorological model. 

10 To ensure the validity of the tropospherically derived data, it is preferred to transmit 
said delay value corrections to a user corresponding to a meteorological temporal 
resolution of said meteorological model information of no greater than 1 hour, and/ or 
corresponding to a meteorological spatial resolution of said meteorological model 
information of no greater than 90 km, insofar as time and distance are linked by speed 

15 of movement of relevant weather features. 

By the above outlined image compression technique it is possible to effect correction 
data transmission to a user at a data rate in the range 25 to 500 bits/s and by 
selectively transmitting only parts of a global image applicable to a user ha a relatively 
small region thereof, to permit correction data transmission effectively at rates well 
2 0 below the top of the range. 

Insofar as the meteorological model derives a tropospheric delay values from data 
employed elsewhere to forecast or predict weather conditions at one or more 
locations, that is, conditions which vary with time, it is possible to predict 
tropospheric delay values in the future from said meteorological information and 
25 develop a prediction set of said corrections for a geographic region of the earth's 
surface, whereby each member of said prediction set describes a correction that 
becomes current as a function of time from development. It is therefore possible to 
communicate said prediction set of corrections as a batch and use members of the set 
as the time for which each was predicted becomes current in respect of the forecast. 
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Such commimication may be to an orbiting satellite and the members re-transmitted 
one at a time as the time for which each was predicted becomes current in respect of 
the forecast. 

A second embodiment employing the method is applicable when the position of the 
5 user receiver with respect to the server and/or GNSS satellites is known. That 
information may be employed by the server with the meteorological information to 
derive actual or mapped tropospheric delay values (rather than zenith delay values) 
for communication to the user for the purpose of setting or correcting the user 
receiver pseudoranging and obtaining accurate timing values. Such communication 

1 0 may be direct or via a network. It may also take place via one or more satellites, such 
as the GNSS satellites as discussed above, although data reduction may be required. 
Insofar as the users location is known, it is not e:xpected to be necessary to derive and 
communicate a set of delay value corrections representing a distribution over a region. 
However, as discussed above, it maybe appropriate to forecast weather conditions for 

15 any user location the user may be in and derive a predicted set of delay corrections 
and communicate these in batch form for use by the user receiver in turn as the time 
for which each member was predicted becomes current. 

According to a second aspect of the present invention there is provided apparatus for 
obtaining data for use by a user of a satelUte positioning system or GNSS, comprising 
2 0 generating means for generating, at a server location remote firom the user from 
meteorological information, at least one accurate tropospheric delay value applicable 
to the user location and means to communicate at least a function of a said value to 
the user as a tropospheric delay correction. 

The server may be arranged to derive a set of tropospheric delay values apphcable to a 
2 5 plurality of user locations. 

In a first embodiment, the apparatus of the preceding paragraph comprises first 
generating means for generating a first set of approximate tropospheric delay values 
firom a first model which is known per se, second generating means for generating a 
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second set of more accurate tropospheric delay values from a said meteorological 
model based on meteorological information, and developing means for developing 
jfrom said first and second delay sets a set of tropospheric delay value modifications 
for use with said first model so that it can provide a set of tropospheric delay values 
5 substantially in agreement with the second set, and said developing means being 
arranged to express the modifications as a set of tropospheric delay corrections. 

Preferably said first generating means utiUses a said first model is based on non- 
meteorological parameters. Also, the developing means may be arranged to express 
said set of corrections each as a diflference betweai corresponding values of the first 
1 0 and second sets, possibly as a firactional change firom the values to be corrected. 

The developing means is arranged to express the corrections as a distribution over a 
region of the earth's surface, preferably in the form of a data file corresponding to a 
greyscale image of multi-bit words, each word representing a location of the region. 
Furthermore the apparatus may include means for compressing said set of corrections. 
15 This may effect lossless compression of the set or for greater reduction, lossy 
compression on the set. 

Each of the first and second generating means may advantageously derive corrections 
for parameters of at least an elevation mapping function used to map the zenith delay 
values to actual delay values. Corrections may be superimposed on the zenith delay 
2 0 correction data set as longer words for communication to the receiver 

The apparatus also transmission means for transmitting said set of corrections to a 
user, and preferably to transmit via an orbiting satellite, which may be a satelUte of 
the GNSS. 

hi a second embodiment, appUcable to apparatus in which the position of a user 
25 receiver of sateUite signals is known, apparatus is arranged to receive firom the user 
information defining at least one of the user location with respect to the server or with 
respect to the GNSS satellites and to provide corrections in the form of tropospheric 
delay values p^ se rather than zenith delay values, although the latter could be 
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provided . 

According to a third aspect of the invention a GNSS user receiver comprises means 
operable to generate from an on-board model from non-meteorological data a set of 
approximate tropospheric delay values q)plicable to identification signals received 
5 from a plurality of said satellites and from and delay values and identification signals 
received from a plurality of said satellites compute an approximate position of the 
receiver relative to the earth's surface or time, means operable to receive a set of 
corrections to said tropospheric delay values derivable from the model, said 
corrections being derived from meteorological data, means to effect modifications to 
10 said derived delay values in accordance with the corrections and means to compute 
the position or time with greater accuracy. 

Said means to effect modification to said delay values may be operable to effect 
interpolation or extrapolation of said corrections according to computed position of 
the user relative to locations for which the corrections have been derived 

15 According to a fourth aspect of the invention a GNSS including a plurality of orbiting 
satellites, apparatus as defined above for obtaining data and a user receiver. 

In the above discussion, tropospheric delay values and zenitihi tropospheric delay 
values have been referred to without regard to their nature. Whereas it is possible to 
derive a single troposphere delay value for a particular position, it is more usual to 
20 derive it as a so-called "wet" delay and a "dry" or "hydrostatic" delay. Apart from 
circumstances where it is important to distinguish, in particular in respect of data 
reduction, in this specification, references to tropospheric delay or delays and their 
derivation is intended to be read as deriving values for each. 

Further details and advantages of the invention will be evident from the following 
2 5 description with reference to the following drawings, in which 

Figure 1 is a schematic representation of a GNSS positioning system known from the 
art, illustrating a xiser positioning receiver device and a plurality of orbiting 
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positioning satellites. 

Figure 2 is a schematic representation of a first embodiment of GNSS positioning 
system embodying the present invention, illustrating a user positioning receiver 
device, a groimd station and a plurality of orbiting positioning satellites, 

5 Figure 3(a) is a graphical illustration of ray-tracing, 

Figure 3(b) is graphical illustration of mapping a reflectivity field in ray tracing. 

Fig 4 is a pictorial representation of a zenith delay data file, suitable for data 
compression by image compression techniques. 

Figure 5 is a graphical representation of how noise affects compression, and 

10 Figure 6 is a schematic illustration, similar to Figure 2 of a second embodiment of the 
invention for a user station at a known position. 

Referring to Figure 1 this represents in schematic form a section of the earths siurface 
50 and in relation thereto a global positioning system 100 comprising a plurality of 
GNSS satelUtes llOi, IIO2, IIO3, ... in earth orbit and a user 120, in the form of a 
15 signal receiver and processor of hand-held or vehicle-mounted type at or above the 
earth's surface. 

The user receiver 120 comprises, in conventional manner, a firont end receiver 122 of 
signals transmitted fi:om the various satellites within view at radio firequency 
frequencies, processing means 124 and information display or like delivery apparatus 
2 0 126. The processing means includes a digital processor that responds to received 
signals fix>m, and characteristic of, the various satellites whose orbital positions are 
known with respect to points on the earth and computes firom the variations in 
reception of these signals a solution comprising the position of the user receiver in 
two- or three-dimensions, and, importantly in some applications, time. 

2 5 The user receiver can determine fix>m the plurality of received satellite signals via the 
above outlined pseudorange measurements an approximate position, but 
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compromised by delays inflicted upon the sig^als by refraction in the troposphere 
caused by refraction, principally the water content of meteorological systems such as 
weather fronts. Such tropospheric refraction may be comp^ated for by, at least to a 
first approximation, by ^plying to processing of the received signals a first model 
130 that represents climatic conditions anticipated for that approximate location at the 
time of year. This so-called climate model is essentially non-meteorological, insofar 
as it is updated infrequently and represents at best a representation of average 
conditions. In known manner the climate model 130 holds parameters for at least one, 
and preferably all of time of year, latitude and altitude as pertain to the position of the 
user and may optionally hold parameters of longitude and/or time of day. 

The climate model is arranged to generate zenith tropospheric delays (ZD) that may 
be applicable to the users location that may be mapped in respect of satellite elevation 
inclination with respect to the receiver to give a more accurate value of tropospheric 
delay and effect pseudorange corrections having regard to the direction actually taken 
15 by the received signal path, particularly if the satellite is at low elevation. 

To this end, the first model may include an elevation mapping ftmction 130' 
employing, for example as a three-term continued fraction approximation, 
substantially as set out by Niell in the paper mentioned above. However, it should 
also be understood that the parameters used in the mapping model, being derived 
2 0 from the time of year, latitude and altitude parameters mentioned above, are also 

subject to errors caused by meteorological disturbances, althou^ for many purposes 
these errors may be considered too small to correct. 

This prior art apparatus utilises the zenith tropospheric delay, and if appropriate 
mapping fimction, to effect an approxinoiate correction to the tropospheric delays 
2 5 enabling the receiver to compute a more accurate solution for position and time. 

Notwithstanding its inh^ent inaccuracies, this first model is valuable insofar as it 
permits a user receiver to be manufactured and used with this (albeit linodted) 
correction &cility built in and not dependant upon receiving signals from elsewhere. 
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Hitherto, the positioning accuracy of a user receiver has been compromised by a 
number of factors but as these resolve, and error sources reduce, it is apparent tiiat 
residual tropospheric delay errors that remain after using the fibrst model are now an 
important cause of limitation to accuracy. 

Referring now to Figure 2, in accordance with the present invention there is provided 
at an earth location a ground station 200, conveniently referred to herein as a server, 
although there may be more than one associated with different regions of the earth's 
surface. This serv^ has no means for receiving satellite signals but is coupled to 
receive from one or more meteorological organisations information representing the 
results of, or suitable for, numerical weather prediction (NWP) for locations at 
various positions around the earth; the meteorological information may be global in 
nature or confined to one or more sub-global regions. 

There is provided in the server, and indicated at 230, a duplicate of the first model 
(130) as used in the user receiver, which contains the aforementioned non- 
1 5 meteorological, climate modelling parameters. 

There is also provided in the server a meteorological model indicated generally at 
250. This model responds to meteorological information provided by the NWP and 
determines accurate values for ZD (as wet and dry components). 

The zenith tropospheric delay ZD values for the two models are compared at 260 in 
2 0 order to determine differences between them that constitutes an error attributable to 
the first model. The differences thus constitute modification values by which the 
product of the first model might be modified or corrected in order to provide the same 
result as the second model. 

These corrections are encoded and subjected to data reduction at 270, as described in 
2 5 detail below, and then communicated to the user receiver by way of a communication 
channel 280 that constitutes an uplink to one or more of the GNSS satellites llOi etc 
by way of transmitter 275 and re-transmission from the satellite or sateUites to the 
user receiver, indicated generally at 220. 
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The user receiver includes a decoder 228 of the correction values data that thus 
provides ZD values for the first model (wet and dry values) effectively making them 
that same as if derived accurately by the second model present only in the remote 
server, for use in the navigation and time computing. 

5 Optionally, as also described below, the corrections may include items applicable to 
mapping functions of the first model so that both zenith delay and msqpping function 
values are given a greater accuracy for the position and time computation. 

The above overview of the system is expanded below with discussion of further 
features that can be employed individually but which when used together interrelate 
10 advantageously. 

Referring to the server 200, the meteorological model relies upon a three-dimensional 
array of grid points for which meteorological information is available and uses such 
information to derive a refiractivity field that permits ray tracing between a ground 
point near the earth's surface and a particular satellite, as a result of which a 
15 tropospheric delay value (for each of the wet and dry delays) can be found. 

At this point it is appropriate to give some background on propagation and 
atmospheric refiractivity and atmospheric effects as they relate to ray-tracing and 
NWP. 

The speed of propagation of an electromagnetic wave through a medium can be 
20 expressed in terms of the refiractive index, «, defined to be the ratio of tiie speed of 
light through free-space to the speed throu^ the medium (Equation 1-1). 

C 

n = - (M) 

V 

where: 

25 n is the refiractive index 

c is the speed of light in fi-ee space 
V is the propagation velocity 
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In practice, and as illustrated in Figure 3(a), the path is curved by refraction as it 
passes between the satellite and earth, most of this delay. The GNSS tropospheric 
time delay, ignoring relativistic effects, is dejSned to be the propagation time of the 
GNSS signal from tihe satellite to the user minus the free space propagation time: 

5 

sv sv 
^trop = ds - \ds . (1^2) 

User User 

where: 

s is the distance along the propagation path. 
10 The first integral is along the curved propagation path; the second integral is 

along a geometric straight path. 



15 



The dififerential equation describing the curved ray path can be expressed, in cartesian 
coordinates, as: 



d ( dr^ 
n 



dsV ds 



= Vn (1.3) 



where r = r(s) is the vector describing the ray path, 
s is the length of the curved ray patii up to r, 
20 n is the refractive index scalar field, 

Vr, a vector field, is the gradient of n. 

The differential equation can be expanded as 

25 

A first order ordinary differential equation (ODE) with known initial values can be 
solved using numerical methods: for example Runge-Kutta or Adams-Moulton 
methods. Higher order differential equations can be solved numerically by rewriting 
30 them as an equivalent system of first order equations. Using the substitution ri = r 
and r2 = r' (the first derivative), the ray path differential equation (1-4) can be 
expressed as an equivalent system of two Gist order differential equations 1-S and 1- 
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6: 




(1-5) 



5 r,' = ^ = ^{Vn-{Vn^:)r:) (1-6) 

as n 

The determination of the ray path therefore amounts to the solution of a system of two 
ODEs with initial values. Standard numerical methods can be used to solve the 
problem: for example, a Runge-Kutta method with adaptive step control consistent 
1 0 with user defined tolerances. 

With the ray path solved, the tropospheric delay can be computed as: 

b c c 

d^ = \n{s)ds+\ds-\ds (i.i) 

a b a 

15 

where a, b and c are as shown in Figure 3(a)* Point b corresponds to the point at 
which ray curvature and refractivity can be assumed to be negligible, in this 
specification above an altitude of 70km. 

2 0 The ray-tracing process to determine the path from user to satellite {a to b in Fig. 3(a)) 
starts at point a and assumes a starting elevation angle of aApparent- Although the 
precise satellite position and therefore aivue is known, aApparent (such that the ray path 
intersects point c) can initially oixly be estimated. Because ray tracing starts off at an 
angle that is at best a guess, the resultant ray path will in general not intersect point c. 

25 By deriving at least two ray traces and using interpolation or iterative methods it is 
possible to establish an angle of suitable accuracy from which tropospheric delay is 
derived. The present invention is predicated upon deriving for use a more accurate 
value for each tropospheric delay. 

30 At least part of tropospheric delay determination is based upon Numerical Weather 
Prediction (NWP) modelling which forecasts the evolution of atmospheric physical 
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processes by applying governing equations, including the conservation of mass, 
momentum and energy. Three-dimensional fields of continuous variables including 
humidity, pressure, temperature and velocity are numerically processed and 
meteorological features, including weather fronts, are secondary derived properties. 
5 A variety of measurements can be input into the numerical model including surface, 
radiosonde and satellite observations. The water cycle is modelled including the 
effects of terrain moisture, sea surface temperature, cloud formation and precipitation. 

Numerical models can be global or of limited area. Limited area high-resolution 
10 models are often termed mesoscale models as they reflect mesoscale meteorological 
features,weather patterns of less than 100km in size. 

The UK Meteorological Office has and makes available so-called Unified Models of 
mesoscale and global data. The NWP model maps each define a grid over the map's 

15 coverage region and the models provide tropospheric delays at corresponding points. 
For example, the UKMO has two NWP models, the so-called global model and the 
mesoscale model. The former has a horizontal resolution of 0.8333 degrees (5/6 
degrees) in longitude and 0.5555 degrees (5/9 degrees) in latitude giving a grid of 432 
X 325 points defining the earths surface, each point associated with a cell of about 

20 60km at mid latitudes and about 90km in the tropics. This global map may be used 
inter alia to provide boundary points for a mesoscale model which is a regional model 
centred on the British Isles and has a resolution of 0.11 degrees in longitude and 
latitude (the grid being rotated with a shiftied pole to maintain uniform horizontal 
resolution) and has 146 x 182 grid points which correspond to an array of cells of 

25 approximately 12km x 12km. Both models have 38 vertical levels and extend to 
about 40km. 

At any grid point of the relevant map The atmospheric refractive index (and therefore 
the gradient of the refractive index) can be derived from numerical weather prediction 
3 0 model pressure, water vapour partial pressure and temperature fields. Atmospheric 
refractivity can be divided into dry (hydrostatic) and wet components. A simple two- 
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term expression having a refractivity accuracy of 0,5% is: 
N-N^^N^^'^-^[p^AUO^ (1-8) 

5 

iV = 10'(/2-l) (1-9) 

where: 

iVis the refractivity. 
n is the refractive index 
10 P is atmospheric pressure (millibar) 

e is the water vapour pressure (milUbar) 
ris temperature (Kelvin) 

The numerical weather prediction fields to be used are expressed in a spherical 
15 coordinate frame, it is computationally convenient, therefore, to generate tiie 
refractivity gradient in spherical coordinates (r,fi,a), which can be converted into local 
curveilinear coordinates (u,v,w) using the following transformation. 

r_ 1 dn \ dn \ dn 

20 L^wJuvw = — ^ + — : Trrv + -— w (i-io) 

-"""^ dr rsmade r da 

A fiuther rotational transformation is then appUed to the local curveiliner coordinate 
frame to give the gradient in a fixed cartesian frame (xj/^z in Fig. 3(b)) suitable for 
numerically solving equations 1-S and 1-6. 

25 

NWP field values between grid-points can be linearly and log-Unearly interpolated. 
Linear and log-linear extrapolation techniques can be used to extend the NWP fields 
beyond the highest grid-point and below the NWP terrain. For an accurate 
constmction of the refiractive index field, account must be taken of the variation of 
30 gravitational acceleration with height and latitude. Hydrostatic equilibrium can be 
assumed. With the three-dimensional refractive index field («) defined, the ray path 
equations can be solved and the tropospheric delay computed using Equation 1-2. 
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The atmospheric refractive index can be divided into dry (hydrostatic) and wet 
components. The wet component possesses the larger spatial and temporal 
variability. It is often convenient to divide total zenith delays into wet and dry zenith 
5 delays. The hydrostatic zenith delay can be accurately modelled given a surface 
pressure measurement, however the wet zenith delay can not be accurately determined 
from surface humidity measurements, as they are not representative of the above 
atmosphere. 

10 In order to aid imderstanding reference is made briefly to discussing the presentation 
of tropospheric delays as so-called "wet" and "dry" components and elevation 
mapping as applied thereto. 

In modelling tropospheric delay, it is convenient in accordance with the first model, to 

relate the tropospheric delay at a given elevation angle, ^, to the zenith delay ( ^Trop ) 
15 by means of a mapping ftmction (m(^)). 

The hydrostatic zenith delay can be accurately modelled given a surface pressure 
20 measurement but the wet zenith delay can not be accurately determined from surface 
humidity measurements, as they are not representative of the above atmosphere. The 
expression for the tropospheric delay at a given elevation angle can be defined as: 

25 

It is noted that care must taken when applying the simplification of the superposition 
of hydrostatic and wet atmospheric delays: the propagation path is dependent upon 
both hydrostatic and wet components. 
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Hydrostatic and wet mapping functions according to the Niell methodology possess a 
high degree of accuracy without the need for prior meteorological information, and 
the variation of tropospheric delay with elevation angle can be efficiently modelled by 
a continued fraction expansion. 

5 

Meteorological features that possess a large spatial and/or temporal variation in 
tropospheric delay will impact the accuracy of NWP-derived tropospheric corrections 
and the bandwidth required for dissemination on a regional or global basis. The 
10 temporal and spatial variation in hydrostatic refractivity is generally small, whereas 
meteorological features associated with rapid changes in atmospheric moisture 
significantly impact the accuracy/bandwidth relationship. 

Meteorological features smaller than the resolution of the numerical prediction model 
1 5 will not be accurately reflected in the NWP-derived tropospheric correction. 

A weather front marks the interface between air masses: defined as a large body of air 
whose physical properties are largely uniform horizontally for hundreds of kilometres. 
The front can mark the occurrence of abrapt changes in atmospheric moisture, 

2 0 temperature and therefore refractivity. Fronts can be divided into three classifications: 

warm, cold and occluded. 

The most rapid change in tropospheric delay is likely to occur when satellite elevation 
and fix)nt inclination are equal. Generally, in the UK, frontal systems move at 30 to 
25 50 kilometres per hour and can result in zenith delay variations of 3 cm/hour. 
Tropospheric delays vary in accordance with inclination to earth insofar as they are 
ray-tracing is taking place through meteorological features that vary differently with 
both altitude and position. 

3 0 The server 200 thus takes as iiq>ut regional or global numerical weather prediction 

model information including pressure, tCTiperature and humidity data and computes the 
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three-dimensional reifractive index field fi-om the meteorological data. Wet and dry 
tropospheric delays are derived for a gridded area, corresponding to the NWP coverage, 
at heights corresponding to a terrain database (that may be the NWP terrain). These are 
transformed to, or initially computed as zenith delays For the same grid locations 
5 including height, wet and dry zenith delays are computed &om the JBrst model 210. 

The server, in computmg the the diflFerence between jSrst model and meteorological 
model, and thus the modifications required to the first model zenith delays to 
makethewm accurate, derives these diferences each as a firactional change firom the first 
1 0 model value, as a percentage. 



There are two benfits to this. Firstly it is found that notwithstanding the actual values of 
the delays and the dififerences, the differences lie in a small range (approximately ± 
10%) 6om flie corresponding first model values; this permits developing a smaller 
15 range of correction values to transmit than if actial value dififerences were used. 
Secondly, it provides for better correction of first model zenith delay values by 
interpolation. 

A correction, -BCj^, defined as a percentage correction to the first model that includes 
20 the variation of zenith delay with height is transmitted. The corrections are in the 
form of a gridded data set. The user can linearly interpolate their tropospheric 
correction firom the adjacent set points. 



BC^(^a.h,) = 100. 

where, 

7T)^ is the zenith delay, measured firom height Ao, computed firom the 

ftwp 

numerical weather prediction, 
30 ZD^^ zenith delay, estimated firom height ho computed using the a 

priori atmospheric model. 



wo 2004/095055 PCT/GB2004/001676 

22 

ho is the height above mean sea level at which the broadcast coxrection 
percentage is computed, 

^is the latitude at which the broadcast coxrection percentage is computed, 
A is the longitude at which the broadcast correction percentage is computed. 

The user can compute the zenith delay as: 



ZD*i,, = ZD^UI^ + BC^yJ 100) 

where, 

BCinm is the interpolated correction percentage, 
h\ is the user height. 



15 It should be noted that this technique avoids requiring the user receiver to store the 
NWP terrain data set, which, as the NWP model evolves, is likely to change. If the 
user is located at one of the grid-points corresponding to the broadcast corrections and 
at hi-hoy the user appUed tropospheric correction is equal to ZD^^ 

20 It will be appreciated that particularly for communicating by satellite and more 
particularly by GNSS satellites there is a limit on transmission bandwidth. Although 
improvements may in the future ease this restriction, for the present it is necessary to 
plan to transmit data at less than 500 bit/s and typically 200 to 250 bits/s. This is 
exacerbated by the satellites not being available to transmit at all times, but only 

2 5 witihdn certain transmission windows. To this end, it is appropriate that the correction 
data derived from the models is suitable for data reduction. 



The NWP model maps each define a grid over the map's coverage region and the 
models provide zenith delays at corresponding points, and thus a set of differences 
3 0 (modifications or corrections) is defined corresponding to said geographical points. 
Thus for a particular set of meteorological information at a particular time, the server 
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can compile a matrix array of such zenith delay modifications for the various 
locations of the organisation's map coverage. 

For example, as mentioned above, the UKMO has two NWP models, the so-called 
global model and the mesoscale model. The former has a resolution giving a grid of 
5 432 X 325 points defining the earths surface, and the mesoscale model which is a 
re^onal model centred on the British Isles has 146 x 182 grid points. 

Thus there exists a two-dimensional matrix array of point correction sets each 
represented by a multi-bit word. In particular this can as in this embodiment be 
represented by an 8-bit word. 

1 0 Thus there exists in the base station, and based upon the particular UM of interest a 
geographical distribution of corrections, essentially an 8-bit greyscale map image of 
the corrections. To the extent that it helps understanding, such a map is capable of 
being represented visually and Figure 4 comprises such a representation of a global 
correction map. 

1 5 Whereas the visual display representation is actually / according to the information 
and its dissemination, it will be appreciated that the format of the information lends 
itself to data compression techniques employed witii such two-dimensional images in 
order to reduce the size of the data file or information for dissemination. 

A lossless or lossy compression map be employed, but in this embodiments the server 
2 0 effects a lossy image compression ; tiie preferred compression is in accordance with 
the JPEG2000 standard (wavelet-based) although other standards, such as JPEG 
(cosine based) or other techniques such as simple sub-sampling may be used, to 
reduce the file size of the infonnation. 

At this point it is appropriate to refer again to the user receiver. Insofar as the 
2 5 correction data signal is received firom the GNSS (or the) satellite with the usual 
signals, no special receiving circuitry is needed. The correction (image) data set is 
passed to processor zenith delay computation which decodes the image for use of the 
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individual pixel values as corrections as described above. Such decoding may be 
accomplished by a hardware feature built into the receiver or such decoding may be 
achieved by software loaded into the central processor of receiver; software for JPEG 
image file decompression is well known 

5 For example, if the file size of the uncompressed image is about 141 kb a 

compression factor of 35 (reducing files size to 4kb) indicates little compression 
noise, but in a greater compression of 140 (to Ikb) compression noise is evident. This 
is also shown graphically in Figure 5. 

It will be appreciated that the conditions for wet delay Amotions are far more complex 
1 0 and produce greater file sizes than dry delay fimctions. Thus separating wet and dry 
delay functions which can reasonably expect compressed file sizes of the order 9 kb 
and 1 kb respectively, there is need to dissenninate about 88,000 bits. 

This data set may be transmitted as a single image to all satellites of a constellation 
group. In this embodiment, the ability to effect data transmission at such a low rate is 
achieved by effecting the transmission to the GNSS satelUtes, with the fiill images at 
high transmission rate / long duration or only part of the image associated with the 
map region associated with a particular satellite or split amongst the satellites for each 
too receive only the part related to it, permitting a fiarther factor of three reduction in 
required reception time; that is, transmission and reception would take about 2 
minutes at a data rate of 250 bits/sec. 

Thus each of the GNSS satellites is able to broadcast, with the normal signals, 
correction signals that each user receiver can employ with the zenith tropospheric 
delay modelling to effect a correction to values employed with the model in 
accordance with substantially current meteorological conditions pertaining to its 
25 location. 

It will be appreciated that the rate of data dissemination is of both technical and 
economic importance. Firstly, the data transmission capacity available for slotting in 
the additional information is limited, at least in current implementations of GNSS. 
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Whereas for other types of data the solution may be to prolong the duration of the 
transmission event, in the case of updating weather dependent data this is not a 
suitable option. Insofar as weather changes and weather features move over the 
earth's surface and thus flie grids employed in the meteorological model, there is a 
5 currency element, that is, a time interval and/or distance for which a desired 

tropospheric delay value is valid. It is believed that such currency time is of tiie order 
of one hour and/or grid size of 50 to 90 km. Thus, is a user is to rely upon 
meteorogically generally tropospheric delay values (by way of a correction) it should 
be within such validity time or position for which generated. Thus, in terms of 
1 0 transmission o a user, transmission must be at such a rate that the user can receive and 
process the information while it is still valid. Satellite transmission rates are both 
slow in bit rate and intermittent in availability for downloading such correction 
information. Thus, there is an imperative to effect a data reduction to accommodate 
transfer of a correction data (image) set whilst the data retains viability. 

1 5 Furthermore, the time taken to download the correction data (image) set should not be 
unduly long as to cause the user to decline to wait for flie time it takes to effect the 
download, decompress the image data file and compute position. 

Thus, it is important to effect the degree of data reduction/image compression that 
2 0 achieves these various objectives. 

Having regard to the above discussion of viability of delay values (and derived 
corrections) it should be noted that employing data based on NWP model approach it 
is in practice possible to predict the weather conditions and derive tropospheric delays 
25 for any given point in advance, up to several hours, notwithstanding that the viability 
of predicted tropospheric delay is relatively small (as described above) once current 
for the predicted time. 

Therefore, it is possible to develop not only a set of tropospheric delays (as Zenith 
30 tropospheric dela>^) associated with a regions grid points but also to develop 
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cotFesponding sets predicted for times in the future. That is, the server cau d^ive a 
prediction set of correction sets. 

This may be useful if the server can otily devote itself periodically to developing delay 
5 corrections, as the prediction set can be stored and its members, correction sets for a 
particular time, be retrieved when that time is current and the set valid. 

Altematively, such a prediction set could be transmitted and stored within the satellite 
for retrieval and re-transmission of correction sets at times for which the validity is 
1 0 current, or analogously stored in tiie receiver for retrieval. 

Notwithstanding the ability to both produce and upload prediction sets of data, in 
practice, at the moment, the data transmission bottle neck is downloading from the 
sateUite (particularly a GNSS satellite) to the user. 

15 

In order to reduce the size of the correction data set filed to be transmitted by a 
satelUte to a user, it may be arranged to transmit only delay corrections applicable to 
users within range, that is, within sight of the satellite, whilst ignoring data for users 
at other points of the globe. This may be achieved by the server determining which 
20 part of the correction data set to upload to any particular satellite, having regard to 
users who can access it, or the server may upload global data set but the satellite 
determine which part of the set to re-transmit. 

The above description has concentrated on describiag the derivation, communication 
25 and usage of corrections to Zenith tropospheric delays which represent major sources 
of positioning errors not correctable by the user receiver first model. 

As mentioned above, the first model also employs a mapping function that employs 
parameters that also depend upon the atmosphere and are accommodated by elements 
30 used in the clinGiate model, that is, a function of time of year (a) latitude (b) and 
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altitude (c). These may also be collected by correction sets derived in the server. This 
moping function may be expressed as: 



5 m(€, a,b,c) = l+ q_ 



1+ b 



l±c (2-1) 

sin(€)+ a 

sin(€)+ b 



10 sin(6)+ c 



where e is the elevation angle. 

15 

This may be further expressed as: 



11 ■ flo+^ 



1, ^>o + A^> 



sm(€) + 



sin(^) + Co + Ac 

where ao, ^^o and cq are the first model values and Aa, A6, Ac are corrections to be 
2 0 applied thereto to effect caapping of elevation derivable by the meteorological model. 



By a fitting process operated such that sum of the squares of the residuals between the 
equation (2-2) and the ray traced (truth) is minimised. Thus a set of correction values 
may be derived composed of Aa, Ab and Ac as a sinoilar data image file or 
25 superimposed upon the Zenith correction data image file by increasing the word 
length thereof to include word elements comprising these corrections. 
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Transmission and reception anti-coding is as before except that the receiver now has 
correction values to employ with the parameters of the mapping function whereby 
both the Zenifli delay and its mat values are more accurately represented. 

5 It will be £^preciated that it is possible to include in elevation mapping functions the 
parameters of longitude and time of day, and the above ^proach to deriving a 
numerical solution therefor may be extended thereto. 

Although it is convenient and in many ways advantageous to employ NWP, other 
10 meteorological soiirces may be used on combination therewith to augment the data 
available. Any other meteorological model that takes a three dimensional view of the 
atmosphere may be employed instead of the NWP. As will be appreciated, the NWP 
derives tropospheric delays are deJBned for a grid having a cell size limited by the 
NWP model in use. Many weather features that have a high moisture content and can 
15 effect tropospheric delay, such a thunder storms, may be below the resolution 
threshold of the NWP model. However, there exists a number of sources of data such 
as satellite images of an essentially dimensional nature that can identify with high 
resolution the existence of such features and the information contained therein can be 
employed to vary the NWP values for a particular cell of the grid to take such features 
2 0 into accoimt across the NWP grid. 

The above described embodiment is intended to enable a user receiver having a built 
in firet, non-meteorological model to determine its position more accurately than is 
possible by use of the model alone. It will be appreciated that part of the computation 
solution is to derive a time values to the accuracy permitted by the model's 
25 mteipretation of tropospheric delay. For some uscts, it is the time function that is of 
importance and such user may know the precise location of a fixed receiver. 

Referring to Figure 6, this shows a schematic representation of a second embodiment 
of GNSS 500. A usCT receiver 520 is similar to the user receiver 120 but lacks (or has 
disabled) the first, non-meteorological model. 
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A server 600 is similar to the server 200 except that the first, non-meteorological 
model is also omitted and transmission is terrestrial rather than via satellite. Within 
the server, tiiere is provided a meteorological processor 650 that uses three- 
dimensional refractive index field generation as described above, based upon NWP 
5 data that may be augmented by additional meteorological data, and effects ray tracing 
capable of deriving tropospheric delay values for a region covered by the NWP data. 
However, the server also receives or has stored thereia, data relating to, or identifying 
the position of the user in earth co-ordinates and satellite elevation and azimuth angle 
with respect to the user. Thus the server does not model a simple zenith delay value 

10 but is able to compute firom the raw data an sqppropriate tropospheric delay value or 
select one of a number of delay values mapped to the correct angles for the user's 
location. The server thus conomunicates the actual delay value to the user by whreless 
or wired means so that the user can derive a more accurate time value and if 
appropriate confirm its position. It will be ^predated that insofar as the 

15 meteorological data is derived for the purpose of forecasting, the derivation of delay 
value and transmission need not be in real time, but could be in advance of its use. 
Not only may a single tropospheric delay value be derived and transmitted this way, 
but in a manner similar to that described for the system 200, a prediction set of delay 
values may be derived and stored, each member of the set being retrieved and used 

2 0 when it becomes current, having regard to the time for which forecast. If the server is 
used for other purposes, fliOT a prediction set may be stored in the server for 
transmission of its members at appropriate times or the set may be transmitted to, and 
stored in^ the user receiver and set members retrieved as tiaey become current. Also, it 
will be ^predated that transmission may be via satellite, using data reduction if 

25 necessary. 

The above disclosure maybe summarised as follows: 



The tropospheric correction server takes as input regional or global numerical weather 
prediction model information including pressure, temperature and humidity data. 
Additional meteorological data, for example high-resolution water v£^our estimates 
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from infrared satellite observations, could be used to augment NWP data and aid in the 
correction of very small meteorological features (for example local convective storms). 

The tropospheric correction server computes the three-dimensional refractive index 
field fix>m the meteorological data. Wet and dry zenith delays are conq>uted for a 
gridded area, corresponding to the NWP coverage, at heights corresponding to a terrain 
database (that may be the NWP terrain). In the case of the server 200, for the same grid 
locations including height, wet and dry zenith delays are computed from a prior model. 
The server then computes the difference between prior and meteorological observation 
based zenith delays as a percentage. The gridded data set of wet and dry percentage 
corrections are then quantised and compressed using a lossy image compression 
technique. The compressed image, including data required for image decoding, is then 
disseminated to the satellite uplink station (part of the satellite navigation system ground 
infrastructure). The image data relates to the current atmospheric state and may also 
include several hours of predicted images. The uploading of data can be from one or 
more uplink stations and can be uplinked to one or more satellites. The data can be 
global or regional in nature. The data sent to each sateUite can be area-limited such that 
the satelUte constellation provides global coverage although each satelUte's data. 

2 0 The batch-uploading to satellites of several hours of predicted images maybe used to 

reduce &e burden on the satellite navigation system ground infrastructure. Each satellite 
broadcasts the most appHcable, most current, tropospheric image as part of the satellite 
navigation system's navigation data. 

25 The user's satellite navigation signal tracking system receives the tropospheric 
images, the standard navigation data including satellite orbit determination parameters 
and makes pseudorange and accumulated carrier observations. The wet and dry 
compressed tropospheric correction images are decompressed. The prior wet and dry 
zenith delays are computed. The prior wet and dry zenith delay model is the same as 

3 0 that used in the tropospheric correction server. The user wet and dry zenith delays. 
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corresponding to the user location, is then computed using interpolation between 
adjacent gridded data points. 

Wet and dry zenith propagation delaj^ are then converted into satellite specific 
5 pseudorange corrections using satellite elevation mapping functions. Finally the 
standard navigation solution is computed, but with tiie addition of meteorological- 
observation-based tropospheric delay corrections. 

A more comprehensive discussion on the factors that affect the implementation of a 
positioning system, numerical weath^ prediction techniques and deriving corrected 
1 0 tropospheric delay information from meteorological measurements, and based may be 
found in the paper " Tropospheric Delay Modelling and Correction Dissemination 
using Numerical Weather Prediction Fields " by M Powe, J Butcher and J Owen 
given in the Proceedings of GNSS 2003 and included in the application from which 
priority for the present application is claimed, and included herein by reference. 
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CLAIMS 

L A method of obtaining tropospheric delay data for use in a satellite positioning 
system or GNSS comprising the steps of generating for a user location, at a 
location remote from the user location and from meteorological information, at 
5 least one accurate tropospheric delay value, applicable to the user location for 
communication as a tropospheric delay correction to a said user. 

2. A method according to claim 1 comprising generating, from a first model which is 
known per se, a first set of approximate tropospheric delay values applicable to 
various user geographical locations, generating from said meteorological model a 

10 second set of tropospheric delay values that are accurate and applicable to said 
various user geographical locations, developing a set of delay value modifications 
for use with said first model so that it can provide a set of tropospheric delay 
values substantially in agreement with the second set, and expressing the set of 
modifications as a set of tropospheric delay corrections for commxmication to a 

15 said user. 

3. A method according to claim 2 wherein ttie first model is based on non- 
meteorological parameters. 

4. A method according to claim 3 wherein said non-meteorological parameters 
comprise at least one of time of year, latitude and altitude. 

20 5. A method according to claim 4 wherein said non-meteorological parameters 
comprise at least one of longitude and time of day. 

6. A method according to any one of claims 2 to 5 wherein said sets of tropospheric 
delay values comprises zenith tropospheric delays. 

7. A method according to claim 6 wherein the first model contains a mapping 
25 fimction relating tropospheric delay at a given elevation angle to the zenith 

tropospheric delay. 
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8. A method according to claim 7 wherein said set of delay value modifications 
comprise a set of modifications for use with the mapping function of the first 
model. 

9. A method according to any one of claims 2 to 8 wherein the modifications are the 
5 differences between corresponding values of the sets attributable to the first and 

meteorological models. 

10. A method according to claim 9 in which the corrections are the modifications. 

11. A method according to claim 9 in which the corrections are the modifications 
expressed as a firactional change from the values of the first set. 

10 12. A method according to any one of claims 2 to 1 1 wherein the set of corrections is 
developed as a data array having values determined for individual grid points on 
the earth's surface and the set of values comprises a distribution of said 
modifications over at least part of the earth's surface. 

13. A method according to claim 12 wherein said set of corrections is derived as a 
15 digital data file. 

14. A method according to claim 13 wherein the set of modifications is derived as a 
set of values Scorn which modifications applicable to points between said grid 
points can be interpolated or extrapolated from said set values. 

15. A method according to claim 13 or claim 14 comprising applying data reduction 
20 to the correction set and deriving a reduced data set for communication to a user. 

16. A method according to claim IS comprising reducing the data size of the 
correction set by an image compression process. 



17. A method according to claim 15 or claim 16 comprising reducing the data size by 
lossy data reduction. 
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18. A method according to claim 17 comprising effecting data reduction by reducing 
the correction set data file according to a JPEG 2000 or JPEG 90 standard. 

19. A method according to claim 1 including providing for said derivation at least one 
of the user position and positions of satellites relative to the user. 

5 20. A method according to claim 19 wherein said derivation at least one of the user 
position and positions of satellites relative to the user is provided by 
communication firom the user. 

21. A method according to any preceding claim wherein said accurate tropospheric 
delay values are derived by a ray tracing technique. 

10 22. A method according to claim 21 wherein said accurate tropospheric delay values 
are derived by three-dimensional refractive index field generation. 

23. A method according to any preceding claim wherein said meteorological model is 
based on nvim©rical weather prediction (NWP) data for a region of the earth. 

24. A method according to any one of claims 21 to 23 wherein said meteorological 
15 information or each said tropospheric delay value correction is augmented by 

directly observed meteorological data. 

25. A method according to claim 24 wherein said directly observed data has a 
resolution smaller than the NWP data. 

26. A method according to claim 24 or claim 25 wherein said directly observed 
20 meteorological data is derived as a data set relating to a region of the earth's 

surface corresponding to at least part of the NWP data. 

27. A metiiod as claimed in any one of claims 23 to 26 wherein said region is 
substantially global. 



28. A method as claimed in any one of claims 21 to 27 comprising predicting 
25 tropospheric delay values in the future from said meteorological information and 
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developing a prediction set of said corrections for said geographic region of the 
earth's siurface, whereby each memb^ of said prediction set describes a correction 
that becomes current as a function of time from developm^t. 

29. A method according to any one of the preceding claims comprising 
5 communicating the correction or set of corrections to a remote user on a 

communication channel or data liiok. 

30. A method as claimed in claims 29 when dependant claim 28 comprising 
communicating said prediction set of corrections as a batch and using members of 
the set as the time for which each was predicted becomes current in respect of the 

10 forecast. 

31. A method according to claims 29 or 30 comprising communicating at least part of 
the corrections to at least one orbiting satellite and re-transmittmg at least part of 
the set to a user from a said orbiting satellite. 

32. A method according to claim 32 comprising communicating the corrections to at 
15 least one orbiting GNSS satellite from which user receives signals to establish at 

least one of position and time. 

33. A method according to claim 31 or claim 32 comprising commxmicating to a said 
sateUite for re-transmission of only that part of the correction data that can be of 
use to a user in a region within range of said satellite. 

20 34. A method according to any one of claims 31 to 33 when dependant on claim 30 
comprising communicating said prediction set of corrections to said at least one 
orbiting satellite in a batch and re-transmitting the members one at a time as the 
time for which each was predicted becomes current in respect of the forecast 



35. A method according to any one of claims 31 to 34 comprising applying data 
25 reduction sufficient to permit transmission of all or part of a said corrections 
useable by a user within a time, dictated by transmission availability and 
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transmission rate of the satellite, substantially lower than the validity time of the 
meteorological information used by the meteorological model. 

36. A melliod according to claim 35 wherein the data reduction is arranged to permit 
transmission of a said corrections to a user corresponding to a meteorological 

5 temporal resolution of said meteorological model information of no greater than 1 
hour. 

37. A method according to claim 35 or claim 36 wherein the data reduction is 
arranged to permit transmission of a said corrections to a user corresponding to a 
meteorological spatial resolution of said meteorological model information of no 

1 0 greater than 50 km. 

38. A method according to any one of claims 35 to 37 wherein the data reduction is 
arranged to permit correction data transmission to a user at a data rate in the range 
25 to 500 bits/s. 

39. A method according to claim 38 wherein the data reduction is arranged to permit 
1 5 correction data transmission in the range 200 to 250 bits/s. 

40. A method of correcting tropospheric delay errors in a GNSS receiver having 
means to derive modelled tropospheric delay values from non-meteorological 
parameters, the method comprising obtaining data in accordance with any one of 
claims 2 to 18 or any one of claims 21 to 39 dependant thereon as a set of 

2 0 corrections, by the use of a first model corresponding to the model in the receiver, 
remotely of the receiver and transmitting said set of tropospheric delay value 
corrections, receiving said corrections in the GNSS receiver and modiJ^ng 
tropospheric delay values obtained by use of the receiver modelling means in 
accordance with the received corrections. 

25 41 . A method according to claim 40 comprising determining within the receiver, from 
the non-meteorological parameters and the model, an approximation to 
tropospheric delay values and therefrom an approximate position of flie receiver 
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relative to the earth's surface, and deriving tropospheric delay value modifications 
from the set of transmitted corrections relevant to the region including said 
approximate position. 

42. A method of detenmning position of a GNSS user receiver having means to 
5 derive modelled tropospheric delay values from non-meteorological parameters, 

comprising determining within the receiver from said non-meteorological 
parameters an approximation to tropospheric delay values and therefrom an 
approximate position of the receiver relative to the earth's surface, deriving 
modified tropospheric delay values in accordance with claim 40 and therefrom 
1 0 determiniag a more accurate receiver position. 

43. A metihod of reducing tropospheric delay errors in a GNSS haviag a user receiver 
that includes means for deriving tropospheric delay values from a non- 
meteorological model, the method comprising deriving, remotely of tiie receiver, 
from a comparable non-meteorological and from a meteorological model a set of 

15 corrections to the values derived from the first model for different locations about 

the earth's stirface, transmitting the set of delay corrections to orbiting satellites 
and re-transmitting the corrections from said satellites to a said user receiver. 

44. A method of providing tropospheric delay values to a receiver at a known 
location comprising deriving at a server location remote from the receiver, using 

20 information relating to the user receiver position and meteorological information, 

at least one tropospheric delay value applicable to the receiver location and 
communicating each said value to the user receiver. 

45. Apparatus for obtaining data for use by a user of a satellite positioning system or 
GNSS, comprising generating means for generating, at a server location remote 

25 from the user from meteorological information, at least one accurate tropospheric 
delay value applicable to the user location and means to communicate at least a 
fimction of a said value to the user as a tropospheric delay correction. 
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46. Apparatus as claimed in claim 45 wherein said meteorological information is 
arranged to derive a set of tropospheric delay values applicable to a plurality of 
user locations. 

47. Apparatus as claimed in claim 46 comprising first generating means for 
5 generating a first set of approximate tropospheric delay values applicable to from 

a first model which is known per se, second generating means for generating a 
second set of more accurate tropospheric delay values from a meteorological 
model based on meteorological information, and developing means for developing 
from said first and second delay sets a set of tropospheric delay value 
10 modifications for use with said first model so that it can provide a set of 
tropospheric delay values substantially in agreement with the second set, and said 
developing means being arranged to express the modifications as a set of 
tropospheric delay corrections. 

48. Apparatus as claimed in claim 47 wherein said first generating means utilises a 
1 5 said fijTst model is based on non-meteorological parameters. 

49. Apparatus as claimed in claim 47 or claim 48 wherein the developing means is 
arranged to express said set of corrections each as a difTerence between 
corresponding values of the first and second sets. 

50. Apparatus as claimed in claim 47 or claim 48 wherein the developing means is 
20 arranged to express said corrections a fractional change from the values to be 

corrected. 

51. Apparatus as claimed in any one of claims 47 to 50 wherein the developing means 
is arranged to express the corrections as a distribution over a region of the earth^s 
surface. 

25 52. Apparatus as claimed in claim 5 1 wherein said region is global. 
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53. Apparatus as claimed in claim 51 or claim 52 in which the developing means is 
arranged to express the corrections in the form of a data file of a greyscale image 
of multi-bit words, each word representing a location of the region. 

54. Apparatus according to any one of claims 47 to 53 including means for 
5 compressing said set of corrections. 

55. Apparatus according to claim 54 in which tihie means for compressing said set of 
corrections comprises means to effect lossy compression on the set 

56. Apparatus according to any one of claims 47 to 55 comprising transmission means 
for transmitting said set of corrections to a user via an orbiting satellite. 

10 57. Apparatus as claimed in claim 56 arranged to transmit to an orbiting satellite 
corrections for re-transmission to the user from that satellite and to transmit only 
corrections applicable to users in regions served by satellite. 

58. Apparatus as claimed in claim 45 or claim 46 wherein the apparatus is arranged to 
receive from the user information defining at least one of the user location with 

1 5 respect to the server or with respect to the GNSS satellites. 

59. A GNSS user receiver comprising means operable to generate from an on-board 
model from non-meteorological data a set of approximate tropospheric delay 
values and fix>m and delay values q)plicable to identification signals received 
from a plurality of said satellites and compute an approximate position of the 

2 0 receiver relative to the earth's surface or time, means operable to receive a set of 
corrections to said tropospheric delay values derivable from the model, said 
corrections being derived from meteorological data, means to ejffect modifications 
to said derived delay values in accordance with the corrections and means to 
compute the position or time with greater accuracy. 

25 60. A GNSS user receiver as clamed in claim 59 wherein said means to effect 
modification to said delay values is operable to effect interpolation or 
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extrapolation of said correctioiis according to domputed position of the user 
relative to locations for which the corrections have been derived 

61. A GNSS including a plurality of orbiting satellites, apparatus as claimed in any 
one of claims 37 to 57 and a user receiver according to claim 59. 

5 62. A GNSS including a plurality of orbiting satellites, apparatus as claimed in any 
one of claims 45 or 58 and a user receiver at a location known with respect to the 
apparatus. 
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